function ind = mapping(x, y, M, N, Xmin, Ymin)
interval = 0.00003;
    ind = (floor((x - Xmin)/interval)* M + floor((y - Ymin)/interval));
end
